Relativistic continuum-continuum coupling in the dissociation of halo nuclei 
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A relativistic coupled-channels theory for the calculation of dissociation cross sections of halo 
nuclei is developed. A comparison with non-relativistic models is done for the dissociation of 8 B 
projectiles. It is shown that neglecting relativistic effects leads to seizable inaccuracies in the 
extraction of the astrophysical S-factor for the proton+beryllium radiative capture reaction. 
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Reactions with radioactive nuclear beams have quickly 
become a major research area in nuclear physics. Among 
the newly developed techniques, the Coulomb dissoci- 
ation method is an important tool to obtain electro- 
magnetic matrix elements between continuum and bound 
states of rare nuclear isotopes pj . These matrix elements 
play an essential role in nuclear astrophysics. At low con- 
tinuum energies, they are the same as the ones involved 
in radiative capture processes of astrophysical interest. 
In particular, the Coulomb dissociation of 8 B projectiles 
allows to extract information on the radiative capture 
reaction p+ 7 Be — > 8 B+7, of relevance for the standard 
solar model and the production of high-energy neutrinos 
in the sun |2(. 

The dissociation of weakly-bound nuclei, or halo nuclei, 
is dominated by the Coulomb interaction, although the 
nuclear interaction with the target cannot be neglected in 
most cases. The final state interaction of the fragments 
with the target leads to important continuum-continuum 
and continuum-bound-state couplings which appreciably 
modify the reaction dynamics. Higher-order couplings 
are more relevant in the dissociation of halo nuclei due 
to their low binding. A known example is the "post- 
acceleration" (or "reacceleration" ]effect observed in the 
dissociation of n Li projectiles 0,0,111 El- 
Two methods have been devised to study higher-order 
effects in projectile dissociation. The method introduced 
in ref. uses the direct solution of the Schrodinger equa- 
tion (DSSE) by space-time discretization. One starts 
with a ground-state wavefunction, propagates it through 
each time-step, and after sufficiently long time the it- 
erated wavefunction is projected into a specific channel 
of interest. Another method discretizes the continuum 
wavefunctions \c) and uses them as input to calculate 
the matrix elements (c'\Vi„t\c} and (c\Vi n t\b), where 
| b) denotes a bound state 0. The matrix elements are 
then used in a coupled-channels calculation for transition 
amplitudes to dissociation channels. This is known as 
continuum discretized coupled-channels (CDCC) method 
and was introduced by Rawitscher [7j to study nuclear 
breakup reactions of the type a + A — > b + c + A. It has 
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been used extensively in the study of breakup reactions 
with stable || and unstable nuclear projectiles 0, U E3 

Special relativity, an obviously important concept in 
physics, is quite often neglected in the afore mentioned 
dynamical calculations. Most rare isotope facilities use 
projectile dissociation at 100 MeV/nucleon. At these en- 
ergies, relativistic effects are expected to be of the order 
of 10%. Relativistic effects are accounted for in the colli- 
sion kinematics, in first-order perturbation calculations, 
but not in dynamical calculations used sofar in the anal- 
ysis of experiments. They also enter the dynamics in a 
non-linear, often unpredictable, way. The reason why 
these effects have not been considered before is due to 
the inherent theoretical difficulty in denning a nuclear 
potential between many-body relativistic systems. Due 
to retardation, an attempt to use a microscopic descrip- 
tion starting from binary collisions of the constituents is 
not possible. A successful approach, known as "Dirac 
phenomen olog y" , has been achieved for nucleon-nucleus 
scattering lllj . But for nucleus-nucleus collisions a rea- 
sonable account of these features has not yet been ac- 
complished. In the present letter, the case of dissocia- 
tion of 8 B projectiles is studied. The major contribution 
comes from the Coulomb interaction with well-known rel- 
ativistic transformation properties. A coupled-channels 
method based on the eikonal approximation with rela- 
tivistic ingredients is developed and compared to semi- 
classical methods 0. Here I omit the consideration of 
other corrections in the eikonal treatment of the scatter- 
ing involving halo nuclei which have also been shown to 
play a relevant role (see, e.g., refs. [f3. IT^I. 

Let us consider the Klein-Gordon (KG) equation with 
a potential Vq which transforms as the time-like com- 
ponent of a four- vector For a system with total 
energy E (including the rest mass M), the KG equa- 
tion can be cast into the form of a Schrodinger equa- 
tion (with h = c = 1), (V 2 + fc 2 - £/) * = 0, where 
k 2 = (E 2 - M 2 ) and U = V (2E - V ). When V < M, 
and E ~ M, one gets U — 2MVq, as in the non- 
relativistic case. The condition \q <C M is met in pe- 
ripheral collisions between nuclei at all collision energies. 
Thus, one can always write U = 2EVq- A further simpli- 
fication is to assume that the center of mass motion of 
the incoming projectile and outgoing fragments is only 
weakly modulated by the potential Vq. To get the dy- 
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namical equations, one discretizes the wavefunction in 
terms of the longitudinal center-of-mass momentum k z , 
using the ansatz 



"J" = S a (z, b) exp (ik a z) 



(€) 



(1) 



In this equation, (z,h) is the projectile's center-of-mass 
coordinate, with b equal to the transverse coordinate. 
4> (£) is the projectile intrinsic wavefunction and (fc, K) 
is the projectile's center-of mass momentum with longitu- 
dinal momentum k and transverse momentum K. There 
are hidden, uncomfortable, assumptions in eq. ^ The 
separation between the center of mass and intrinsic coor- 
dinates is not permissible under strict relativistic treat- 
ments. For high energy collisions we can at best justify 
eq. nfor the scattering of light projectiles on heavy tar- 
gets. Eq. \\\is only reasonable if the projectile and target 
closely maintain their integrity during the collision, as in 
the case of very peripheral collisions. 

Neglecting the internal structure means 4>k a (£) = 1 
and the sum in eq. ^rcduces to a single term with a — 0, 
the projectile remaining in its ground-state. It is straight- 
forward to show that inserting eq. ^ m the KG equa- 
tion (V 2 + k 2 - 2EVu) * = 0, and neglecting V 2 5 (z, b) 
relative to ikdzSo (z, b), one gets ikdzSo{z 1 h) = 
EVoSo (z,h), which leads to the center of mass scatter- 
ing solution <So(z,b) = exp —iv~ 1 f^_ oo dz' Vo(z',h) , 

with v = k/E. Using this result in the Lippmann- 
Schwinger equation, one gets the familiar result for 
the eikonal elastic scattering amplitude, i.e. fo = 
—i (k/2n) J db exp (zQ • b) {exp [ixO 3 )] — 1} > where the 
eikonal phase is given by exp[z%(b)] = <So(oo,b), and 
Q = K' — K is the transverse momentum transfer. 
Therefore, the elastic scattering amplitude in the eikonal 
approximation has the same form as that derived from 
the Schrodingcr equation in the non-relativistic case. 

For inelastic collisions we insert eq. \I\m the KG equa- 
tion and use the orthogonality of the intrinsic wavefunc- 
tions (j>k a (£). This leads to a set of coupled-channels 
equations for S a : 



k 2 ) S a e lk ° 



with the notation la) 



(2) 



Neglecting terms of the 



form V 2 S a (z, b) relative to ikdzS a (z, b), eq. EJreduces 
to 

. v dS a (z,b) = ^ {a]Vola/) Sa/{Zjh) e i( fcQ ,-M*. (3 ) 



The scattering amplitude for the transition — > a is 
given by 

/«(Q) = -^ J dh exp(zQ-b) [S a (b) - S at o] , (4) 



with S a (b) — S a (z — oo, b). The set of equations [3] and 
Elare the relativistic-CDCC equations (RCDCC). 

I have used the RCDCC equations to study the dis- 
sociation of 8 B projectiles at high energies. The ener- 
gies transferred to the projectile are small, so that the 
wavefunctions can be treated non-relativistically in the 
projectile frame of reference. In this frame the wave- 
functions will be described in spherical coordinates, i.e. 
| a) = \jlJM), where j, I, J and M denote the angular 
momentum numbers characterizing the projectile state. 
Eq. [3] is Lorentz invariant if the potential Vq transforms 
as the time-like component of a four-vector. The matrix 
element (a \Vq\ a 1 ) is also Lorentz invariant, and we can 
therefore calculate them in the projectile frame. 

The longitudinal wavenumber k a ~ (E 2 — M 2 ) 1 / 2 also 
defines how much energy is gone into projectile exci- 
tation, since for small energy and momentum transfers 
k' a — k a = (E' a — E a ) /v. In this limit, eqs. and 0] re- 
duce to semiclassical coupled-channels equations, if one 
uses z = vt for a projectile moving along a straight- 
line classical trajectory, and changing to the notation 
S a (z,b) = a a (t,b), where a a {t,b) is the time-dependent 
excitation amplitude for a collision wit impact parameter 
b (see eqs. 41 and 76 of ref. ^3). Here I use the full 
version of eq. 0] 

If the state |a) is in the continuum (positive 
proton+ 7 Be energy) the wavefunction is discretized ac- 
cording to \a;E a ) = JdE' a T(E' a ) \a;E' a ), where the 
functions T(E a ) are assumed to be strongly peaked 
around the energy E a with width AE. For convenience 
the histogram set (eq. 3.6 of ref. Q) is chosen. The 
inelastic cross section is obtained by solving the RCDCC 
equations and using da / dfldE a = \f a (Q)|~ T 2 (E a ). 

The potential Vq contains contributions from the nu- 
clear and the Coulomb interaction. The nuclear po- 
tentials are constructed along traditional lines of non- 
relativistic theory. The standard double-folding approx- 
imation F^- aT ^(R) = J p a (r) vq (s) pt (r') is used, where 
vq (s) is the effective nucleon-nucleon interaction, with 
s = r + R r'. The ground-state densities for the pro- 
ton, 7 Be (p a ), and Pb target (pr), are taken from ref. 
[Tlj . The M3Y effective interaction 0] is used for 
vq (s). The nucleus- nucleus potential is expanded into 
1 = 0, 1, 2 multipolarities. These potentials are then 
transformed as the time-like component of a four-vector, 
as described above (see also ref. ^3)- The multipole 
expansion of the Coulomb interaction in the projectile 
frame including retardation and assuming a straight-line 
motion has been derived in ref. |14| . The first term 
(monopole) of the expansion is the retarded Lienard- 
Wiechert potential which does not contribute to the ex- 
citation, but to the center of mass scattering. Due to 
its long range, it is hopeless to solve eq. with the 
Coulomb monopole potential, as S a (z, b) will always 
diverge. This can be rectified by using the regulariza- 
tion scheme S a (b) — > exp [i (2?yln (kb) + \ a {b))\ S a (b), 
where ij = ZpZx/hv and the S a (b) on the right-hand 
side is now calculated without inclusion of the Coulomb 



3 



"O 



CD 
"D 

To 

"D 

OJ 



1.00 = 



0.10 



0.01 



0.100 r ''/ 



0.010 




2 4 6 8 

6 [degrees] 



FIG. 1: Angular distributions for the dissociation reaction 
8 B+Pb -> p+ 7 Be+Pb at 50 MeV/nucleon. Data are from 
ref. [2^ and are corrected for the detection efficiency e. The 
dotted curve is the first-order perturbation result of ref. [25|. 
The solid curve is the RCDCC calculation. The dashed curve 
is obtained with the replacement of 7 by unity in the nuclear 
and Coulomb potentials. 



monopole potential in eq. [3J The purely imaginary ab- 
sorption phase, Xa (b), was introduced to account for ab- 
sorption at small impact parameters. It has been cal- 
culated using the imaginary part of the "tpp" interac- 
tion ^t|, with the 8 B density calculated with a modified 
Hartree-Fock model ^8|. The El and E2 interactions, 
taken from ref. 0] replacing vt — > z. Explicitly, 



El/j. 



{b 2 + 7 2 z 2 ) 3/2 
for the El (electric dipole) field, and 



T-& if M = ±1 
V2z if fx = , 

(5) 



/ 3"7T 

v E2fl = ^/-eY 2fl u 



b 2 if p = ±2 
T-27 2 6z if p = ±1 
273 (2 7 2 z 2 - b 2 ) if 



(6) 



H = 



for the E2 (electric quadrupole) field, where e\ = |e and 
e 2 = §§e are effective charges for p+ 7 Be. For 7 — > 1 
these potentials reduce to the non-relativistic multipole 
fields (see, e.g., eq. 2 of ref. [ijj) in distant collisions. 

A single-particle model was used for 8 B with a Woods- 
Saxon potential adjusted to reproduce the binding energy 
of 0.139 MeV 11 H2 . I follow the method of ref. H 
and divide the continuum into bins of widths AE a = 100 
keV, centered at E a = 0.01, 0.11, 0.21, 1.01 MeV, bins 
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FIG. 2: Cross sections for the dissociation reaction B+Pb — > 
p+ 7 Be+Pb at 83 MeV/nucleon and for 9s < 1.8°. Data are 
from ref. [2^ ]. The dotted curve is the first-order perturbation 
result. The solid curve is the RCDCC calculation. The dashed 
curve is obtained with the replacement of 7 by unity in the 
nuclear and Coulomb potentials. 



of widths AE a = 250 keV, centered at E a = 1.25, 1.5..., 
2.0 MeV, and bins of width AE a = 0.75 MeV, centered 
at E a = 2.50, 3.25, 10.00 MeV. Each state a is a 
combination of energy and angular momentum quantum 
numbers a = {E ai l, j, J, M}. Continuum s, p, d and / 
waves in 8 B were included. 

The calculations with the RCDCC equations were com- 
pared to the data of refs. [24J and [23. The results were 
folded with the efficiency matrix as well as the energy 
averaging procedures explained in ref. [24| and provided 
by the RIKEN collaboration 0. At 83 MeV/nucleon, 
the angular distribution was chosen to match the same 
scattering angles referred to in ref. [23l |. 

Figure ^ shows the angular distributions for the 
dissociation reaction 8 B+Pb— +p+ 7 Be+Pb at 50 
MeV/nucleon. Data arc from ref. The relative 

energy E between the proton and 7 Be is averaged over 
the energy intervals E=0.5-0.7 MeV (upper panel) and 
E = 1.25 — 1.5 MeV (lower panel). The dotted curve is 
the first-order perturbation result reported in ref. [25| . 
The solid curve is the result of the RCDCC calculation. 
The dashed curve is obtained with the replacement of 7 
(Lorentz factor) by unity in the nuclear and Coulomb 
potentials. The dashed curve is on average 3-6% lower 
than the solid curves in figure 1. Using non-relativistic 
potentials yields results always smaller than the full 
RCDCC calculation. It is a non-trivial task to predict 
what the relativistic corrections do in a coupled-channels 
calculation. 

Figure 2 shows the relative energy spectrum between 
the proton and the 7 Be after the breakup of 8 B on lead 
targets at 83 MeV/nucleon. The data are from ref. [23| . 
In this case, the calculation was restricted to b > 30 fm. 
The dotted curve is the first-order perturbation calcu- 
lation, the solid curve is the RCDCC calculation, and 
the dashed curve is obtained with the replacement of 7 
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by unity in the nuclear and Coulomb potentials. The 
difference between the solid and the dashed-curve is of 
the order of 4-9%. I have also used the DSSE method, 
described in ref. [(|, to compute the same spectrum, 
with the same partial waves, assuming a large lower cut- 
off in the parameter of b — 30 fm. This would justify 
a semiclassical limit. The same relativistic nuclear and 
Coulomb potentials have been used in both calculations. 
The difference between the two results (the RCDCC and 
the DSSE) is very small (less than 2%) for the whole 
range of the spectrum. To my knowledge such a compar- 
ison has never been made before. This is a proof that the 
two methods yield the same result, if the same potentials 
are used, and as long as large b's are considered. Such a 
comparison is only possible because the same potential 
was used for the bound-state and the continuum. The 
DSSE method |(| does not allow for use of different po- 
tentials for p+ 7 B. This is not a problem in the RCDCC 
method, since the states \b) and |c) can be calculated 
within any level of sophistication, beyond the simple po- 
tential model adopted here. In this respect, the RCDCC 
is superior than the DSSE method and is more suitable 
for an accurate description of dynamical calculations. 

The conclusions drawn in this work are crucial in 
the analysis of Coulomb breakup experiments at high 
bombarding energies, as the GSI experiment at 254 
MeV/nucleon [26j. In table 1 I show the calculations for 
the correction factor A — i^RCDCC 



T CDCC 



CDCC 



for the dissociation of 8 B on lead targets at 3 bombard- 
ing energies. E is the relative energy of the proton and 
7 Be. One sees that the relativistic corrections tend to in- 



crease the cross sections. At 250 MeV/nucleon they can 
reach a 15% value. This has been treated before in first- 
order perturbation theory, but not in the dynamical cal- 
culations with continuum-continuum coupling used in the 
experimental analysis 23,26]. The consequence of using 
these approximations on the extracted values of the as- 
trophysical S-factors for the 7 Be(p,7) reaction in the sun 
is not easy to access. It might be necessary to review the 
results of some of these data, using a proper treatment 
of the relativistic corrections in the theory calculations 
used in the experimental analysis. Other improvements 
of the present formalism needs to be assessed. The rela- 
tivistic effects in the nuclear interaction has to be studied 
in more depth. The effect of close Coulomb fields H^IH^ 
should also be considered in the case of dissociation of 
halo nuclei. 



Lab energy 


A 


A 


A 


[MeV/nucleon] 


E = 0.1 MeV 


E = 1 MeV 


E = 2 MeV 


50 


1.5% 


4.2% 


3.4% 


80 


3% 


5.5% 


4.1% 


250 


5.3% 


14.6% 


6.9% 



Table 1: Relativistic corrections in the dissociation of 
8 B projectiles impinging on lead targets at different bom- 
barding energies. E is the relative energy of the proton 
and 7 Be. 
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